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@ Single-Point Access to Data Distributed on Many Processors 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A description of the functions and data 
structures is defined that would be neces- 
sary to implement the Chapel concept of 
distributions, domains, allocation, access, 
and interfaces to the compiler for trans- 
formations from Chapel source to their 
run-time implementation for these con- 
cepts. A complete set of object-oriented 
operators is defined that enables one to 
access elements of a distributed array 
through regular aridnnetic index sets, giv- 


ing the programmer the illusion that all 
the elements are collocated on a single 
processor. This means that arbitrary re- 
gions of the arrays can be fragmented and 
distributed across multiple processors with 
a single point of access. This is important 
because it can significantly improve pro- 
grammer productivity by allowing the pro- 
grammers to concentrate on the high- 
level details of the algorithm without 
worrying about the efficiency and commu- 


nication details of the underlying repre- 
sentation. 

This work was done by Mark James of Caltech 
for NASA’s Jet Propulsion Laboratory. Further in- 
formation is contained in a TSP ( see page 1 ). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-42505. 


© Estimating Dust and Water Ice Content of the Martian Atmosphere 
From THEMIS Data 


NASA’s Jet Propulsion Laboratory, Pasadena, California 


Researchers at JPL and Arizona State 
University conducted a comparative 
study of three candidate algorithms for 
estimating components of the Martian 
atmosphere, using raw (uncalibrated) 
data collected by the Thermal Emission 
Imaging System (THEMIS) . THEMIS is 
an instrument onboard the Mars 
Odyssey spacecraft that acquires image 
data in five visible and nine infrared 
(IR) wavelength bands. The algorithms 
under study used data collected from 
eight of the nine IR bands to estimate 
the dust and water ice content of the 
atmosphere. Such an algorithm could 
be used in onboard data processing to 
trigger other algorithms that search for 
features of scientific interest and to re- 
duce the volume of data transmitted to 
Earth. 


The algorithms studied were based 
on regression models. In the study, 
the optical depths estimated by these 
algorithms were compared with opti- 
cal depths estimated in ground-based 
processing using fully calibrated data 
from both THEMIS and the Thermal 
Emission Spectrometer (TES). TES is 
an instrument onboard the Mars 
Global Surveyor spacecraft that also 
observes the planet at infrared wave- 
lengths, but at a lower spatial resolu- 
tion than THEMIS does. Of the algo- 
rithms studied, the one that 
performed best was based on a Gauss- 
ian Support Vector Machine regres- 
sion model. The test results indicated 
that this algorithm, operating on the 
raw data, had error rates that were 
within the uncertainty associated with 


the estimates obtained by the ground- 
based analysis of the fully calibrated 
data. This level of fidelity demon- 
strates that these algorithms are suffi- 
ciently accurate for use in an onboard 
setting. 

This work was done by Kiri Wagstaff Re- 
becca Castaho, and Steve Chien of Caltech 
for NASA’s Jet Propulsion Laboratory and 
Joshua Bandfield of the Arizona State Uni- 
versity. Further information is contained in 
a TSP (see page 1 ). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPQ-43590. 


© Computing a Stability Spectrum by Use of the HHT 

Unlike in the predecessor method, the mathematical sign of the damping is retained. 

Goddard Space Flight Center, Greenbelt, Maryland 


The Hilbert-Huang transform (HHT) 
is part of the mathematical basis of a 
method of calculating a stability spec- 
trum. This method can be regarded as 
an extended and improved version of a 
prior HHT-based method of calculating 


a damping spectrum. In the prior 
method, information on positive damp- 
ing (which leads to stability) and nega- 
tive damping (which leads to instability) 
becomes mixed into a single squared 
damping loss factor. Hence, there is no 


way to distinguish between stability and 
instability by examining a damping spec- 
trum. In contrast, in the present stability- 
spectrum method, information on the 
mathematical sign of the damping is re- 
tained, making it possible to identify re- 
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gions of instability in a spectrum. This 
method is expected to be especially use- 
ful for analyzing vibration-test data for 
the purpose of predicting vibrational in- 
stabilities in structures (e.g., flutter in air- 
plane wings) . 

A brief summary of the HHT is pre- 
requisite to a meaningful brief summary 
of the present method. The HHT has 
been a topic of several prior NASA Tech 
Briefs’ articles, the first and most detailed 
being “Analyzing Time Series Using 
EMD and Hilbert Spectra” (GSC-13817), 
NASA Tech Briefs, Vol. 24, No. 10 (Octo- 
ber 2000), page 63. To recapitulate: The 
HHT method is especially suitable for 
analyzing time-series data that represent 
nonstationary and nonlinear physical 
phenomena. The method involves the 
empirical mode decomposition (EMD), 
in which a complicated signal is decom- 
posed into a finite number of functions, 
called “intrinsic mode functions” 
(IMFs) , that admit well-behaved Hilbert 
transforms. The HHT consists of the 
combination of EMD and Hilbert spec- 
tral analysis. 

An unavoidably lengthy description of 
the mathematical basis of the prior 
damping-spectrum method is also pre- 
requisite to a meaningful brief summary 
of the present method. The instanta- 
neous amplitude of a vibration signal at 
time t is given by 


n 

x(t)='£ l Cj(t) + r n 
7=1 

where n is an integer, cf t) is an IMF, and 
r„ is a residue signal. 

For each IMF (for example, the kth 
one), a Hilbert transform is performed 


to obtain a complex time-dependent 
function: z k (t) = c k (t) + id h ( t ) . 

The time-dependent amplitude 
[a h {t)}, phase [0*(t)], and frequency 
[®j(t)] of the Ath IMF are then given by 
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a k (t) = [ck(t) + d£(t)\ ; 

6 k (t )= tan -1 and 
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Q} k (t) = - 


dt 


The damping of the Mi IMF is given 
by 


. , 2 da k (t) 
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a k ( t ) dt 

The damping loss factor of the Mi IMF 
is then given by 
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Then summing all the squared damping 
loss factors as functions of time and fre- 
quency and letting frequency become a 
continuous variable CO, one obtains the 
damping spectrum rf(co,t), which is re- 
lated to an amplitude spectrum a(co,t) 
via the equation 

-2 da(co,t) 
a(co,t) dt 


T]\co,t) = - 


-.2 


® 2 + 


1 d.a(co,t) 
a{(0,t) dt 


This concludes the description of the 
prior method. 

In the present method, one computes 
a damping loss factor rh,(t) or r; ( co,t ) by 
use of equations similar to those shown 
above, but with the following notable 
differences: 

• Instead of using the Hilbert transform 
to compute a complex function and 
then using the complex function to 
compute the amplitude function, one 
uses a cubic spline to compute the am- 
plitude function. The reason for this 
change is that in a practical implemen- 
tation, a Hilbert transform can intro- 
duce spurious oscillations that can mask 
true damping or anti-damping, whereas 
any spurious oscillations introduced by 
a cubic spline are much smaller. 

•The instantaneous frequency CO k (t) or 
t o(t ) is not calculated as indicated above. 
Instead, it is calculated by use of the nor- 
malized HHT. This change is necessitated 
by a limitation of the Hilbert transform — 
too complex to discuss here — that has 
been a topic of prior publications. 

• One retains the sign of the damping by 
simply refraining from squaring the 
damping loss factor: in other words, 
r\ ( CO,t ) becomes the stability spectrum. 
Areas of positive and negative damp- 
ing can be readily distinguished on a 
plot of the spectrum. To make areas of 
negative damping even more readily 
apparent, it could be desirable, in 
some cases, to place areas of positive 
damping and areas of negative damp- 
ing on separate plots. 

This work was done by Norden E. Huang 
of Goddard Space Flight Center. Further in- 
formation is contained in a TSP ( see base 1 ). 
GSC-14833-1 


26 


NASA Tech Briefs, March 2007 


